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SUMMARY 

Infection with dengue, the most prevalent mosquito-borne virus, manifests as dengue fever (DF) 
or the more fatal dengue haemorrhagic fever (DHF). DHF occurs mainly when an individual 
who has acquired antibodies to one serotype is inoculated with another serotype. It was reported 
that mosquito control may have increased the incidence of DF and DHF due to age-dependency 
in manifesting these illnesses or an immunological mechanism. Tetravalent dengue vaccine is 
currently being tested in clinical trials. However, seroconversions to all four serotypes were 
achieved only after three doses. Therefore, vaccines may predispose vaccinees to the risk of 
developing DHF in future infections. This study employed an individual-based computer 
simulation, to emulate mosquito control and vaccination, incorporating seroconversion rates 
reported from actual clinical trials. It was found that mosquito control alone would have 
increased incidence of DF and DHF in areas of high mosquito density. A vaccination 
programme with very high coverage, even with a vaccine of suboptimal seroconversion rates, 
attenuated possible surges in the incidence of DF and DHF which would have been caused 
by insufficient reduction in mosquito abundance. DHF cases attributable to vaccine-derived 
enhancement were fewer than DHF cases prevented by a vaccine with considerably high 
(although not perfect) seroconversion rates. These predictions may justify vaccination 
programmes, at least in areas of high mosquito abundance. In such areas, mosquito control 
programmes should be conducted only after the vaccination programme with a high coverage 
has been initiated. 

Key words: Dengue haemorrhagic fever, mathematical modelling, vaccine safety, vaccines, 
vector control. 
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in Singapore [6]. As the average age of dengue infec- 
tion is increasing in many countries [7], this mechan- 
ism, known as 'endemic stability' [8], may explain 
the rising incidence of DF in some countries. 

The unique immunological aetiology of DHF may 
have led to a similar unwanted consequence as well 
[9]. DHF is known to occur in a secondary infection 
more frequently than in a primary infection, due 
to immunological mechanisms including antibody- 
dependent enhancement (ADE) [10]. On the other 
hand, many of the individuals in areas of high vector 
mosquito abundance would be infected by, and acquire 
immunity against, multiple serotypes while they are 
clinically protected by this cross-immunity [11]. 
Consequently, these individuals develop resistance to 
DHF unknowingly, since those infected by more than 
one serotype rarely manifest DHF [12]. As mosquito 
abundance decreases, an increasing number of individ- 
uals would experience secondary infections after the 
protective cross-immunity has waned, and the incidence 
of DHF would increase. Consistent with this hypo- 
thesis, incidence of DHF was in a negative relationship 
with mosquito abundance in Thailand [9, 13]. 

Whether these seemingly paradoxical hypotheses 
are true or not, the concern that insufficient mosquito 
control may increase the incidence of DF and DHF 
underscores the importance of understanding the 
effect of mosquito control and future vaccination on 
the target population. In effect, two tetravalent live 
vaccines were developed and tested in clinical trials: 
a classical live-attenuated vaccine from the Walter 
Reed Army Institute of Research (WRAIR) [14] and 
a chimeric yellow fever-dengue composite vaccine 
(ChimeriVax or CYD) from Sanofi-Pasteur [15]. Of 
these, only CYD remained in phase II b trials [16]. 
However, dengue vaccines harbour a unique concern: 
enhancing antibodies may be induced by the vaccine 
itself [17, 18]. In particular, a classical live-attenuated 
vaccine, which induces immune responses similar to 
those induced by wild-type viruses [19], may predis- 
pose a vaccinee to the risk of enhancement upon a 
subsequent infection. The risk of inducing enhance- 
ment by CYD, which is a yellow fever vaccine frame- 
work inserted with prM and E genes from dengue 
viruses, is regarded as being less than for classical 
live-attenuated vaccine [19]. However, 'only phase 
IV trials and post-marketing surveillance will provide 
a definitive answer as to whether ADE constitutes a 
risk for vaccinees' [15]. Therefore, seroconversion to 
all four serotypes is regarded as a prerequisite for a 
tetravalent dengue vaccine. By contrast, in clinical 



trials conducted to date, live tetravalent vaccines 
induced antibodies to less than four serotypes in a con- 
siderably large proportion of vaccinees, even after two 
consecutive injections [20, 21]. Although serocon- 
versions to all four serotypes were achieved after 
three injections, such a three-dose regimen of a live 
vaccine is unprecedented. Furthermore, administering 
all three doses to all vaccinees may be difficult to 
achieve in developing countries. 

It is recognized that mathematical models are 
useful to predict the population-level effects of dengue 
vaccine [22]. Alternatively, the present study used 
an individual-based model, based upon results from 
clinical trials. Mathematical models creates a set of 
differential equations. In contrast, individual-based 
models create a large number of human individuals 
in the computer's memory, and observes their behav- 
iour [23]. Diverse scenarios, regarding seroconversion 
rates of dengue vaccine and vaccination coverage 
were compared. In addition, the incidence of DHF 
attributable to ADE derived from prior vaccination 
was estimated. Finally, the optimal strategy for den- 
gue control is discussed. 



MATERIALS AND METHODS 

Assumptions for protective and enhancing antibodies 

While the titre of antibodies against dengue virus 
is high, the protective role of antibodies is dominant: 
however, as the titre wanes, antibodies enhance de- 
velopment of DHF [24]. It is not known whether pro- 
tective antibodies and enhancing antibodies are 
physically separable (Fig. la) ([25] and E. Konishi 
et ai, personal communication), or whether the 
same antibodies switch from protection to enhance- 
ment as their titres wane — in other words, virtual 
enhancing antibody (Fig. lb) [26]. However, these 
two hypotheses converge into an identical software 
coding if the first hypothesis (Fig. la) surmises that 
the enhancing antibody persists life-long and if the 
second hypothesis (Fig. lb) assumes that antibody 
induced by viral inoculation exerts life-long enhance- 
ment. This framework enabled the distinction between 
enhancing antibodies induced by prior wild-type infec- 
tion and those attributable to vaccine. 

An inoculation with a wild-type virus induces pro- 
tective antibodies specific to that serotype. On the 
other hand, an individual inoculated with a tetra- 
valent live vaccine acquires protective antibodies specific 
only to the serotype(s) to which s/he seroconverted. 
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Fig. 1. Two hypotheses regarding antibodies to dengue 
viruses, (a) Enhancing antibodies and protective antibodies 
are different. Enhancing antibodies react with all virus 
serotypes (broken arrows). In contrast, protective antibodies 
are specific to a serotype (Dl in this figure), but exert 
transient cross-reactive protection against other serotypes 
(solid arrows), (b) The same antibodies may play different 
roles in protection and enhancement, depending on their 
titres. Antibodies (specific to Dl in this figure) exert pro- 
tection against other serotypes when their titres are high 
(solid arrows). As the titre wanes, these antibodies act as 
enhancing antibodies (broken arrows). Dl, D2, D3, and D4 
represent dengue virus serotypes 1, 2, 3, and 4, respectively. 

The protective antibodies are assumed to exert life- 
long serotype-specific protection, as well as transient 
cross-serotype protection for a duration of C years 
from the latest inoculation. Sabin [11] observed volun- 
teers who were inoculated sequentially with two differ- 
ent strains of dengue virus, with the interval between 
the inoculations being maximally 9 months. It was 
found that cross-protection against severe illness still 
persisted at least after this interval. In addition, vacci- 
nation with yellow fever-dengue 2 chimeric vaccine 



induced cross-serotype protection that lasted >1 year 
[27]. Furthermore, it was reported that secondary 
infections resulted in DHF, DF, and asymptomatic 
infections at 2-6, 1-9, and 1-6 years after the primary 
infections, respectively [28]. These observations sug- 
gested that cross-serotype protection against DHF 
may last for >1 year. Therefore, C was assumed to 
follow a normally distributed probability distri- 
bution function (PDF) with a mean equal to 2 years 
(Table 1). An individual inoculated with either wild- 
type or vaccine acquires enhancing antibodies, which 
persist throughout that individual's life [29, 30]. 

Assumptions for DF and DHF 

Individual-based model simulation software (detailed 
in Protocol SI of reference [13]) was modified to 
describe immunological behaviour of the host 
(Fig. 2a). When a naive individual is inoculated by 
wild-type dengue virus, s/he transitions to the cross- 
protected state. In the course of this transition, s/he 
may develop DF with an age-dependent probability 
defined in Figure 2b (constructed based upon refer- 
ence [5] and P. G. Coleman, personal communi- 
cation), and DHF with a fixed small probability of 
0-2% [12]. If an individual in the cross-protected 
state is inoculated with a virus serotype, s/he acquires 
antibodies specific to this serotype and remains in the 
cross-protected state. C years after the most recent 
inoculation, the individual moves to the expired cross- 
protection state. When an individual in this state is 
inoculated with a serotype to which s/he does not pos- 
sess specific antibodies (i.e. unexperienced serotype), 
s/he may manifest DF. The individual may also 
develop DHF with an enhanced probability of 4% 
[12], if s/he already possesses enhancing antibodies. 
An individual who has seroconverted to 'L' serotype 
transitions to the completely immune state (Table 1). 
The transmissibility (or viraemia) may be enhanced 
'T'-fold during manifestation of DHF [31] (Table 1). 

Basic and effective reproductive numbers 

In the mathematical models for dengue proposed 
so far, transmission intensity was often expressed as 
mosquito density [32, 33] or as basic reproductive 
number (R 0 ) [9, 34, 35]. In the present study, R 0 rep- 
resents transmission intensity, since R 0 is proportional 
to vector abundance [3, 36, 37]. The range of R 0 was 
selected, considering previous estimates [34, 38, 39]. 

The immunological state of each individual was 
updated at discrete time-steps of 2 weeks' length, 



1248 U. Thavara and others 



Table 1 . Variable parameters given to simulations 



Definition 


Symbol 


Unit 


Values assumed 
in simulations 


Cross-protective period 


c 


Year 


v YJr (mean — z, s.d. — u d) 


im umoer 01 serotypes necessdry 10 comer 




Lyimensioniess 


ivdnuoniiy seiecieu irom z 9 j or *f 


complete resistance to DHF 








Enhancement of transmissibility during 


T 


Dimensionless 


PDF (mean =10, s.d. = 2-5) 


manifestation of DHF 








Basic reproductive number 


Ro 


Dimensionless 


Fixed* (range 1-20) 


Inhomogeneous mixing 


I 


Dimensionless 


Fixed* (range 0-1) 


Vaccination coverage 


V 


% 


Fixed* (range 0-100) 


Seasonality 


S 


Dimensionless 


Fixed* (range 0-1) 


Total fertility rate 


TFR 


Dimensionless 


Fixed* (range 0-1) 



PDF, Probability distribution function with normal distribution; s.d., standard deviation. 

* Fixed, a fixed value was input into each simulation to examine the effect of the parameter on the simulation result. 



which approximates the sum of intrinsic incubation 
and infectious periods [11, 31, 40]. The risk of being 
infected by a serotype (force of infection or viral 
inoculation rate) at the zth time-step (Fi) was obtained 
as: 

Fi = R 0 ,i-i x Ui-x/Ni-u (1) 

where Rqj-i, K-i and N f _i represent basic repro- 
ductive number, viraemic load and population size 
in the (/— l)th time-step, respectively. Here, the virae- 
mic load at the zth time step (U t ) is defined as: 

Ui = DF, + (T x DHF,), (2) 

where DF, and DHF, denote the number of patients 
with DF and DHF, respectively. T is the above- 
mentioned enhancement in transmissibility (Table 1). 
Effective reproductive number, which is defined as 
the secondary infectious cases originating from a 
primary infectious case, was estimated for each time 
step to be compared with R 0 given to the simulation. 

Estimation of vaccine seroconversion rates 

To exclude the effect of natural infections, the results 
of clinical trials conducted in non-endemic areas were 
used in the present study (Table 2). Two trials adopted 
a regimen of 0-3-5-12 months [20, 21]. However, the 
short first interval resulted in trans-serotype inter- 
ference [15, 21]. Therefore, the present study assumed 
0-6-12 months interval (Table 3). The results in adult 
vaccinees were applied to paediatric vaccinations. All 
the trials in Table 2 defined a 50% plaque reduction 
neutralization titre (PRNT 50 ) of 1:10 as the cut-off 
for seropositivity. This titre had been considered as 
the surrogate for maximal protection, as in Japanese 
encephalitis [41]. The clinical studies reported the 



seroconversion rate after the first dose (RmvX that 
after the second dose CRm+2%) and that after the 
third dose (^,1 + 2+3%) for each serotype i (7=1, 2, 3, 
or 4). From these values, the seroconversion rate 
achieved solely by the second dose (^,2) was esti- 
mated as: 

R U2 = 100 xCRf.i+2 -^,i)/(100-^,i), (3) 

Similarly, the seroconversion rate achieved solely by 
the third dose (^,3) was estimated as: 

R h3 = 100 x 1+2+3 -i?M+ 2 )/(100 - RiA+ll (4) 

Here, it was assumed that vaccinees that had serocon- 
verted would not turn seronegative. 

Compromised protective seroconversion rates 

Rij (i= 1—4, 7=1—3) represents the maximal achiev- 
able protective potencies. However, the efficacies 
reported in the clinical trial were low (i.e. 30-2%) [16], 
indicating that 'protective' seroconversion rates for 
CYD may be much lower than R ij9 which was 
measured based on a PRNT 50 of 1:10. Therefore, the 
present study conducted a sensitivity analysis by con- 
sidering tetravalent vaccines with compromised protec- 
tive potencies (Table 3). CYD 40 and CYD 2 o were 
assumed to induce protective seroconversion rates of 
40% and 20%, respectively, compared to the serocon- 
version rates measured as PRNT 50 of 1:10. CYD 10 o 
induces protective seroconversion rates equal to those 
measured as PRNT 50 of 1:10. 

Vaccine assumptions 

Enhancing antibodies derived from vaccine and 
those from wild-type virus were differentiated in the 
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Fig. 2. Individual-based model for dengue infections, (a) Diagram of the transition between immunological states. 
Transition between immunological states (filled arrow) occurs as a result of either viral inoculation or expiration of time 
from the most recent inoculation (open arrow). During a state transition, the individual may manifest dengue fever (DF) 
or dengue haemorrhagic fever (DHF) (arrow head). * A modification of the model in reference [13] was made: the 
individual is predisposed to the risk of DHF only if enhancing antibodies pre-exist, (b) Age-dependent probability for 
an infection to manifest as DF in an individual who is not protected specifically or cross-reactively. The probability is 
expressed as: 100/[1 + l/exp(— 3-44 + 0-177 x age)] % ([5] and P. G. Coleman, personal communication). 



software. Inoculation by a vaccine will not manifest as 
clinical illness or lead to development of transmissible 
viraemia. Since CYD tetravalent vaccine contains 
only two dengue-virus peptides (E and prM), this 
vaccine may not confer efficient cross-protection. 
Consistent with this assumption, the efficacy of 
CYD measured within 2 years after vaccination was 
only 30% [16]. Therefore, this study adopted an un- 
favourable scenario for CYD, assuming that CYD 
does not induce cross-serotype protection. 



Emulation of vaccination programme 

A future vaccination programme was emulated. At the 
100th year in a 150-year simulation, a vaccination 
programme was initiated. After this year, an attempt 
was made to vaccinate all children who reached 
age 2 years but only with a successful coverage of 
V%. This age was selected because it was the youngest 
age tested in the clinical trials so far conducted. It was 
attempted to re-vaccinate the vaccinees at 6 months 
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Table 2. Serotype-specific seroconversion rates (%) of CYD tetravalent 
dengue vaccine based on a PRNT 50 of 1:10, obtained from clinical trials 
conducted in non-endemic areas 





Age 


Time of 
vaccination 


Dl 


D2 


D3 


D4 


Ref. [20], Figure 4a 


2-45 












After 1st dose 




0 month 


15 


32 


50 


67 


After 2nd dose 




3-5 months 


50 


75 


72 


87 


After 3rd dose 




12 months 


77 


85 


90 


89 


Ref. [21], Figure 2a 


18-45 












After 1st dose 




0 month 


12 


67 


27 


64 


After 2nd dose 




4 months 


70 


93 


73 


87 


After 3rd dose 




12-15 months 


100 


100 


100 


100 


Ref. [21], Figure 2b 


18-45 












After 1st dose 




0 month 


30 


72 


33 


88 


After 2nd dose 




8-11 months 


92 


100 


100 


100 



Dl, D2, D3, and D4 represent dengue virus serotypes 1, 2, 3, and 4, respectively. 



Table 3. Dengue serotype-specific seroconversion rates 
(%) assumed in simulations 





Dl 


D2 


D3 


D4 


Seroconversion rates from [21] Figure 2a 








30 


72 


33 


88 


Ri,\+2 


92 


99 


99 


99 


Ri,\+2+?> 


100 


100 


100 


100 


CYD 100 










*U 


30 


72 


33 


88 


Ri,2 


89 


96 


99 


92 


^,3 


100 


100 


100 


100 


CYD 40 










*u 


12 


29 


13 


35 




36 


38 


40 


37 


^,3 


40 


40 


40 


40 


CYD 20 












6 


14 


7 


18 




18 


19 


20 


18 


^,3 


20 


20 


20 


20 



and i?f,i + 2+3 denote the seroconversion rate for 
serotype i after the first, second, and third doses, respect- 
ively. R i2 and R i 3 represent the seroconversion rate achieved 
solely by the second and third doses, respectively. R i2 and 
R i3 of CYD 10 o were estimated using equations (3) and (4) 
in the main text, respectively. R^, R i2 , and R i 3 of CYD 40 
and those of CYD 2 o were 40% and 20% of the correspond- 
ing values for CYD 10 o, respectively. R i2 and R i3 were 
input into the simulations. Dl, D2, D3, and D4 represent 
dengue virus serotypes 1, 2, 3, and 4, respectively. 



and 12 months after the first dose, but with a success- 
ful follow-up coverage of V%, respectively. This was 



because 100% accomplishment of any vaccination 
programme, especially of three-dose regimen, is diffi- 
cult to achieve in developing countries. 

Inhomogeneous mixing 

The degree of inhomogeneity in the mixing pattern 
was represented by /, where /= 0 indicated completely 
homogeneous mixing and /= 1 indicated maximally 
inhomogeneous mixing (see Supplementary online 
Appendix 1 and Table 1). 

Seasonality 

Since vector mosquito abundance fluctuates season- 
ally [42], seasonality is an important factor which 
should be considered in modelling a mosquito-borne 
disease [33]. Therefore, the present study expressed 
seasonality as a sinusoidal function of R 0 , as follows: 

R 0 (0) = R 0 x [1 - cos(6>) x S], (5) 

where 0 represents the phase of a time-step in a 
year [0^#^27r], Rq{0) indicates R 0 at this phase, 
and Ro expresses the average annual R 0 . S represents 
the strength of seasonality. S= 1 indicates the strong- 
est seasonality and S = 0 represents no seasonality 
(Table 1). 

Age-specific birth rate and mortality rate 

Demographic data from Thailand were used as an 
example of a region heavily afflicted by dengue. 
Age-population structure of Thailand in 1960 was 
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Fig. 3. Results from simulations plotted over mosquito abundance and vaccination coverage. Mosquito abundance 
is represented as basic reproductive number (R 0 ). Vaccination coverage is defined as V in the Methods section, (a-c) 
Viral inoculation rate (/1000 individuals per year), incidence (/100000 individuals per year) of dengue fever (DF) (d-f), 
dengue haemorrhagic (DHF) (g-i), and DHF attributable to vaccine-derived antibody-dependent enhancement (j-l), 
were averaged from the last 30 years in each 150-year simulation and then averaged from 20 simulations, (a, d, g, j) 
CYD 10 o, {b, e, h, k) CYD 40 , and {c, f, i, I) CYD 2 o are compared. Parameter setting for the simulations was (TFR, I, S) = 
(2, 0, 0-2). 



1252 U. Thavara and others 



(a)R 0 = 4 




0 20 40 60 80 100 



(b)R 0 = ^5 

400 | 




o I J 

0 20 40 60 80 100 





+ CYD 100 CYD 40 --- --- CYD, 



Fig. 4. Incidence of dengue fever (DF), dengue haemorrhagic (DHF), and DHF attributable to vaccine-derived antibody- 
dependent enhancement (ADE), from simulations that assumed moderate or high mosquito abundance, plotted over 
vaccination coverage. Incidence (/l 00 000 individuals per year) of (a, b) DF, (c, d) DHF, and (e, f) DHF attributable to vaccine- 
derived ADE, were compared in CYD 10 o, CYD 40 , and CYD 2 o- Simulations were conducted assuming (a, c, e) moderate 
mosquito abundance (R 0 = 4) or (b, d, f) high mosquito abundance (Ro = 15). Parameters were set to the same values as in 
Figure 3. 



used as the initial population structure. Population 
growth was represented by the total fertility rate 
(TFR) [43]. Based on the TFR given to a simulation, 
age-specific birth rate was reconstructed for the each 
age class of females (Supplementary Appendix 2). The 
age-specific mortality rate reported from Thailand 



in 2005 was used throughout the entire simulation. 
TFR of Thailand was 1-86 in 2000. 

Source code 

The source code, which was written in PERL language, 
can be obtained from Supplementary Appendix 3. 
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Fig. 5. Effects of inhomogeneous mixing (7), seasonality (S), and total fertility rate (TFR) at different levels of vaccination 
coverage. Results from simulations are plotted over vaccination coverage and (a, d, g) inhomogeneous mixing (/), (b, e, h) 
seasonality (S), and (c, f, i) TFR. R 0 was set at 15. (a-c) Viral inoculation rate (/1000 individuals per year), incidence 
(/l 00 000 individuals per year) of (d-f) dengue fever (DF) and (g-i) dengue haemorrhagic (DHF) were averaged from the 
last 30 years in each 150-year simulation and then averaged from 20 simulations. CYD 40 was used. Parameter settings 
were: (TFR, S) = (2, 0-2) for (a, d, g); (TFR, I) = (2, 0) for (b, e, h); (I S) = (0, 0-2) for (c, / /)• 



RESULTS 

Effects of vaccination and vector control on the 
incidence of DF and DHF 

As a result of simulations, the viral inoculation rate 
(or force of infection) decreased linearly as the mos- 
quito abundance (represented by R 0 ) decreased or 
as the vaccination coverage increased (Fig. 3 a-c). 
However, the behaviour of the incidence of DF and 
DHF was more complex. A steep high ridge of DF 
incidence occurred (Fig. 3d-f), while the ridge of 
DHF incidence was blunt (Fig. 3 g-i). DHF incidence 
decreased substantially as the vaccination coverage 
increased at any given mosquito abundance 



represented by R 0 (Fig. 3 g-i). However, reduction 
in mosquito abundance would not necessarily 
decrease the incidences unanimously. In particular, 
at low vaccine coverage (<60%), reducing the mos- 
quito abundance from a high level (R 0 > 15) to a mod- 
erate level (^0*^4) would increase the incidence, 
especially of DF. The maximal incidence of DHF 
attributable to vaccine-derived ADE was in the 
order of CYD 20 , CYD 40 , and CYD 100 (Fig. 3j-l). 
This incidence of vaccine-derived DHF for CYDioo 
was relatively small at a very high vaccination cover- 
age (Fig. 3j). In contrast, the incidence of vaccine- 
derived DHF for CYD 40 and CYD 2 o was highest at 
100% vaccination coverage (Fig. 3k, I). The effect of 
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Fig. 6. Temporal patterns in simulations employing different dengue control strategies. Incidence of dengue fever (DF) and 
dengue haemorrhagic (DHF) is plotted over years in simulation. The strategy was based on either (a, d, g) vector mosquito 
control only (b, e, h, j, I) vaccination only, or {c, f, i, k, m) vaccination followed by vector control. The temporal change in 
mosquito abundance (R 0 ) and vaccination coverage are presented in panels (a-c), while incidence (/l 00 000 individuals 
per year) is presented in (d-f, j, k) for DF and in (g-i, /, m) for DHF. CYD 40 was used in {e, f, h, i) while CYD 2 o was used in 
(/-/, m). Parameters were set to the same values as in Figure 3. 



vaccination on effective reproductive number is pre- 
sented in Supplementary Appendix 4. 

Preventive and predisposing effects of vaccines 

Figure 3 was intersected at R 0 = 4 (moderate mosquito 
abundance) and i?o = 15 (high mosquito abundance; 



Fig. 4). The incidence of DF responded to vaccination 
coverage differently between high and moderate mos- 
quito abundances: DF incidence was more refractory 
to vaccination in areas of high mosquito abundance 
(Fig. 4b) than in areas of moderate mosquito abun- 
dance (Fig. 4a). The vaccines affected DHF incidence 
in a contrasting manner. At both high and moderate 
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mosquito abundance, all the vaccines reduced 
DHF incidence (Fig. 4c, d). Vaccine-derived DHF 
(Fig. 4e,f) was much less frequent than DHF pre- 
ventable by CYD 10 o and CYD 40 (Fig. 4c, d), while 
this margin was more obscure for CYD 2 o- 

Effect of inhomogeneous mixing, seasonality, and 
population growth 

The effects of inhomogeneous mixing, seasonality, 
and population growth rate on the epidemiological 
parameters were examined at different levels of vacci- 
nation coverage (Fig. 5). It is intuitive that less 
inhomogeneous mixing and larger population growth 
were associated with increased viral inoculation rate 
(Fig. 5a, c). Interestingly, stronger seasonality was 
correlated with slightly increased viral inoculation 
rate (Fig. 5b). However, the response of disease inci- 
dence to increased viral inoculation rate was counter- 
intuitive. Under the situation of high mosquito 
abundance (R 0 =\5), the increased viral inoculation 
rate led to decreased DF incidence (Fig. 5 d-f), as 
predicted by the endemic stability hypothesis. The 
incidence of DHF and its relationship to vaccination 
coverage was affected only slightly by the variation 
in inhomogeneous mixing, seasonality, and popu- 
lation growth rate (Fig. 5 g-i). 

Temporal patterns 

In Figure 6, the temporal patterns in the incidence 
of DF and DHF are presented for three control strat- 
egies: mosquito vector control only (Fig. 6a), vacci- 
nation only (Fig. 6b), and vaccination followed by 
vector control (Fig. 6c). In the area of moderate mos- 
quito abundance, all three strategies reduced the inci- 
dence of both DF and DHF (data not shown). In 
contrast, in an area of high mosquito abundance 
(R 0 =15), mosquito control of moderate achievement 
(to Ro = 4) led to an increase in incidence of DF 
(Fig. 6d) and DHF (Fig. 6g). Vaccination alone by 
CYD 40 did not affect DF incidence markedly 
(Fig. 6e), but reduced DHF incidence noticeably 
(Fig. 6h). Vaccination using CYD 40 followed by vec- 
tor control achieved a substantial reduction in inci- 
dence of both DF (Fig. 6f) and DHF (Fig. 6i). 
Vaccination alone by CYD 2 o did not exert a notice- 
able influence on the incidence of DF (Fig. 6j) or 
DHF (Fig. 61). However, vaccination by CYD 2 o 
which preceded vector control attenuated the potential 
increase in incidence that would have resulted from 




Fig. 7. Optimal strategy to reduce viral inoculation rate and 
incidence of dengue fever (DF) and dengue haemorrhagic 
(DHF). The optimal strategy, which reduces the viral 
inoculation rate and incidence of DF and DHF, is 
superimposed on the results of simulations which used 
CYD40 (Fig. 3). This strategy, which is represented by the 
curved arrow in (a) and (b) is composed of a vaccination 
phase, and a mosquito control phase. Initially, by 
attaining high coverage of vaccination in 'vaccination 
phase', (a) DF incidence and (b) DHF incidence decrease. 
In the subsequent 'mosquito control phase', R 0 is reduced, 
thereby decreasing these incidences to a lower level. With 
this strategy, combining vaccination and mosquito control, 
the ridges of incidence of (a) DF and (b) DHF can be 
circumvented. 

vector control alone for DF (cf. Fig. 6d,k) and DHF 
(cf. Fig. 6g,m). 

DISCUSSION 

Individual-based model simulation has become 
increasingly common for comparing disease control 
strategies [44^6]. The present study employed this 
methodology to predict the effects of vector control 
and vaccination on the incidence of dengue-related 
diseases. 
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As a result, it was predicted that mosquito control 
alone is likely to increase the incidence of both DF 
and DHF in areas of high mosquito abundance. 
Although extremely strong suppression of vector 
abundance could decrease the incidence eventually, 
incidence would increase transiently while mosquito 
reduction remains incomplete. Despite this concern, 
mosquito control activities will probably be continued 
in developing countries, partly because being endemic 
for dengue virus reduces the attractiveness for tourism 
and overseas investment. Since the vector mosquito 
larvae infest intended water containers as well as 
disposals [47, 48], not only mosquito control but 
also improvement in water supply and garbage col- 
lection systems will decrease mosquito abundance 
[49, 50]. 

The present study emulated seroconversion rates of 
CYD, the dengue vaccine with the highest prospect of 
proceeding to phase III trial. Sensitivity analyses were 
conducted to compensate for the uncertainties in sero- 
conversion rates. Consequently, the main conclusion 
was not affected by these uncertainties; the vaccine 
could attenuate the possible surge in DF and DHF 
driven by the decrease in mosquito abundance. 
However, the predisposing effect of vaccine-derived 
ADE would vary greatly depending on protective 
seroconversion rates. It should be noted that the least- 
favourable assumptions for a vaccine were adopted, 
including that vaccines always induce enhancing anti- 
bodies, the enhancing capability persists throughout 
life, and the vaccine does not confer cross-protection. 
Hence, the estimates of a predisposing effect of vac- 
cine may be exaggerated. In addition, DHF cases 
caused by vaccine-derived ADE would be much 
fewer than DHF cases prevented by a vaccine of con- 
siderably high (although not perfect) seroconversion 
rates. Therefore, a vaccination programme which 
uses such a sub-optimal vaccine may be justifiable, 
at least in areas of high mosquito abundance 
where mosquito reduction may increase the incidence 
of DF and DHF. As represented by the arrows in 
Figure 7, the peaks of DF/DHF incidence could be 
circumvented if mosquito control is preceded by 
high vaccination coverage. Since some may feel that 
any vaccine that predisposes vaccinees to the risk of 
DHF cannot be ethically acceptable, the present 
study provides quantitative information to ethical 
and economic discussions of this issue. Although vac- 
cination of not only small children but a large part 
of the population may be necessary at the initial 
phase of a vaccination programme [51], the present 



study did not investigate this important topic. 
Further studies are warranted on these issues. 

Collectively, the present study proposes a new 
methodology to predict and compare the population- 
level effect of dengue vaccines. The prediction can be 
updated easily as seroconversion rates are improved, 
or as currently unknown parameters are reported 
from field/experimental studies. The predictions 
made here, however peculiar they may appear, should 
be considered in developing a global dengue control 
strategy. 
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